Nonlinear dynamos at infinite magnetic Prandtl number 
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The dynamo instability is investigated in the limit of infinite magnetic Prandtl number. In this 
limit the fluid is assumed to be very viscous so that the inertial terms can be neglected and the flow 
is slaved to the forcing. The forcing consist of an external forcing function that drives the dynamo 
flow and the resulting Lorentz force caused by the back reaction of the magnetic field. The flows 
under investigation are the Archontis flow, and the ABC flow forced at two different scales. The 
investigation covers roughly three orders of magnitude of the magnetic Reynolds number above onset. 
All flows show a weak increase of the averaged magnetic energy as the magnetic Reynolds number is 
increased. Most of the magnetic energy is concentrated in flat elongated structures that produce a 
Lorentz force with small solenoidal projection so that the resulting magnetic field configuration was 
almost force-free. Although the examined system has zero kinetic Reynolds number at sufficiently 
large magnetic Reynolds number the structures are unstable to small scale fluctuations that result 
in a chaotic temporal behavior. 

PACS numbers: 05.45.-a, 91.25.Cw 



Dynamo instability is considered to be the main mech- 
anism for the generation and sustainment of magnetic 
field throughout the universe [l[ . In this scenario an ini- 
tially small magnetic field is amplified by currents in- 
duced solely by the motion of a conducting fluid. This 
process saturates when the magnetic field becomes strong 
enough for the Lorentz force to act back on the flow and 
reduce its ability to further amplify the magnetic energy. 
The exact effect the Lorentz force has on the flow and 
which property of the flow is altered to prevent further 
magnetic field amplification is the subject of many cur- 
rent investigations. 

One can in principle envision many scenarios for sat- 
uration especially in the presence of turbulence. For ex- 
ample, saturation can occur because the magnetic field 
suppresses the fluid flow thus reducing the magnetic 
Reynolds number of the flow to its critical value. This 
behavior is expected close to the dynamo onset for lam- 
inar flows. In another scenario, dynamo can saturate by 
suppressing the chaotic stretching of the magnetic field 
lines 0, Q ; or the dynamo can saturate because the fold- 
ing of the field lines is modified to be less constructive. 
Finally, saturation of the dynamo can be due to correla- 
tions between the velocity and the magnetic field without 
suppressing the ability of the flow to amplify an uncor- 
rclatcd passive vector field thus in principle still be a 
dynamo flow (see for example [H-Q)- 

In principle there is no unique answer to this ques- 
tion and the list above does not exhaust all possibili- 
ties. The exact mechanism can be a combination of dif- 
ferent processes whose choice can depend on the type 
of forcing and the examined parameter range. Identify- 
ing the mechanisms in simple flows can help in unveiling 
the dependence of the amplitude of magnetic field sat- 
uration on the magnetic Reynolds number R m and the 
fluid Reynolds number Re. This dependence is impor- 
tant for astrophysics, bearing in mind that both num- 



bers are very large, with 1 <IC Rm <C Re for the Sun and 
1 <C Re <C Rm for the galaxy. 

In this large parameter space, it is interesting to exam- 
ine special limits of the governing equations, in order to 
obtain a basic understanding of the non-linear dynamo 
behavior. This is attempted in this work. Here satu- 
rated dynamos are investigated in the simplified limit of 
large viscosity. In this system the inertial terms can be 
neglected and the resulting flow is a Stokes flow driven 
by the external force and the Lorentz force. This limit 
is typically used as a model for the small scale galactic 
dynamo for which the magnetic Prandtl number (which 
is equal to the ratio of the viscosity to the magnetic dif- 
fusivity) is very large. In previous investigations [rj-TToj 
the forcing was assumed random and changing in time in 
order to model turbulent fluctuations at the viscous cut- 
off scale. The turbulent fluctuations at this scale result 
in the fastest growing dynamo modes at the kinematic 
stage of the dynamo. Here the opposite limit is examined 
and a time independent forcing is used. A steady forcing 
might not be suitable in describing the turbulent velocity 
fluctuations, however it allows to investigate the reaction 
of the magnetic field to large scale flows with long time 
correlations that could be more relevant in the saturated 
stage. Another interesting property of this limit is that 
the system is stripped off all hydrodynamic instabilities. 
Thus, when comparing with results of turbulent MHD 
flows one can distinguish which effect is due to turbulence 
and hydrodynamic instabilities and which solely due to 
the effect of the Lorentz force. 

This paper is structured as follows, the next section de- 
scribes the dynamical equations in the limit considered. 
Section II explains the numerical method and its limi- 
tations. Section III describes the results for each of the 
examined flows. Section IV gives some common scaling 
laws observed for all flows, and conclusions are drawn in 
the last section. 
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I. EQUATIONS AND SET UP 

We consider flows in a triple periodic domain of size 
L = 2ir driven by a simple time independent body force. 
The non-dimensional MHD equations for a unit density 
fluid are then given by 

Pe(<9 f u + u • Vu) = b Vb- VP + V 2 u + F (1) 
dtb + u- Vb = b Vu+ — V 2 b. (2) 

Rm 

F is the external body force that in the absence of mag- 
netic fields sustains a velocity field u, of unity root mean 
squared value U = (u 2 ) 1 / 2 = 1.0 (where (•) stands for 
spatial average). The Reynolds number is defined as 
Re = UL/{2i:u) where v is the kinematic viscosity and 
we are going to consider the limit Re -C 1. The mag- 
netic field is expressed by b = \/Reb/U, where b is 
the dimensional magnetic field measured in units of ve- 
locity. The factor of \[Re has been introduced to set 
the Lorentz force (that leads to saturation) to the same 
order of magnitude as the forcing term. The magnetic 
Reynolds number is defined as R m = UL/(2wq), where r/ 
is the magnetic diffusivity. P is the pressure that ensures 
the incompressibility condition V • u = 0. 

We are interested in the limit Re <C 1 while R m is still 
finite; thus the magnetic Prandtl number P m , defined as 
the ratio of the two Reynolds numbers P m = Re/R m = 
tends to infinity. In this limit then, the left hand 
side of equation Q] can be dropped and velocity field is 
slaved to the forcing and the Lorentz force and is given 
by 

u = -V~ 2 [b • Vb - VP + F] , (3) 

where V~ 2 stands for the inverse laplacian. Eq. |3] repre- 
sents the Stokes- flow of a viscous fluid under the influence 
of the body-force F and the solenoidal projection of the 
Lorentz force b • Vb — VP. The approximations made 
here, are valid in the small Reynolds number limit pro- 
vided that the velocity gradients do not become too sharp 
or the velocity time scale too short. In the examined 
calculations the velocity field appears to have a smooth 
behavior both in time and space as R m is increased, so 
this condition is not violated. 

Equations [2] and [3] form a closed set of equations that 
are under investigation in this work. They are going to 
be referred to as the infinite magnetic Prandtl number 
equations. The evolution equation for the magnetic en- 
ergy is then given by: 

k(6 2 ) = (b.(Vu).b)--^((Vfc) 2 ) (4) 

Multiplying cq. [3]by V 2 u space averaging and subtract- 
ing from eq. [4] we obtain 

U(b 2 ) = (u • F) - ((Vu) 2 ) - ^<(V6) 2 ). (5) 



Eq. [5] is the energy balance equation that expresses that 
the rate of change of magnetic energy is equal to the rate 
energy is injected to the system by (u • F) minus the 
energy lost by the viscous damping e u = ((Vu) 2 ) and 
Ohmic dissipation e& = ((S/b) 2 ) / R m . It worth pointing 
out that even in the limit R m — > oo energy is dissipated 
at all scales due to the viscous term. It is also noted 
here that the kinetic energy of the flow does not enter 
the energy balance equation since u is a slaved vector 
field. If we did add it, it would be of order Re smaller 
than the magnetic energy. In the remaining text we will 
refer to the kinetic and magnetic energy as E u = ^(u 2 ) 
and Eb = |(b 2 ) respectably, although for a true "dimen- 
sional" comparison between the two a factor of Re should 
be included (ie E tot = E u + E b /Re). 

The infinite magnetic Prandtl equations conserve one 
ideal invariant (for R m = oo): the magnetic helicity fill ]. 
Multiplying equation [2] by the vector potential a (such 
that b = V x a) and space averaging we obtain 

k(a.b) = -l(b.Vxb>, (6) 

^ rim 

the magnetic helicity H m = | (a ■ b) is thus conserved for 
infinite R m . For large but finite R m this invariance is ex- 
pected to be violated at scales i ~ 1/ \/R m . Since there 
is no large scale source term for magnetic helicity, it can 
be generated or destroyed only at these small scales. In 
dynamo investigations the initial magnetic field is very 
small and thus also the initial magnetic helicity is negli- 
gible. However a flow with positive (negative) hydrody- 
namic helicity will generate negative (positive) magnetic 
helicity in the larger scales where it will "pile-up" and 
positive (negative) helicity in the small scales where it 
will be dissipated. Thus although the flow itself does not 
generate helicity it drives preferentially one sign of helic- 
ity to the small scales where it is dissipated. So if the 
saturated state is magnetically helical with one sign of 
helicity it implies that the opposite sign of helicity has 
been dissipated at the small scales. 

II. NUMERICAL METHOD 

The infinite magnetic Prandtl number equations 12 131 
were solved in a triple periodic domain of size L = 2tt us- 
ing a standard pseudo-spectral method and a third order 
in time Runge-Kuta. The code was based on a full MHD 
code [HI, Il3| The truncated fields were dialliased every 
time a quadratic nonlinear term was calculated based on 
the 2/3 rule. 

The resolution used varied from 64 3 up to 256 3 grid 
points. With a resolution of 256 3 one can typically 
achieve magnetic Reynolds numbers of the order of a few 
hundreds for the kinematic problem depending on the 
flow. However, in the non-linear regime, the velocity am- 
plitude is significantly reduced (in some cases by two or- 
ders of magnitude) and the magnetic energy is dissipated 
also by the viscous term making the flow less efficient at 



3 



generating small scales. Thus, although at the examined 
resolution only magnetic Reynolds numbers of few hun- 
dreds can be achieved for the kinematic problem, for the 
nonlinear problem magnetic Reynolds numbers of a few 
thousand can easily be achieved if the initial conditions 
were also at the nonlinear stage. The following procedure 
was used for the numerical runs that were performed: 
Starting from a well resolved kinematic dynamo simula- 
tion (very small initial magnetic field) the system was 
evolved until non-linear saturation was achieved. The 
output of this run was then used as initial conditions for 
the next run at a higher magnetic Reynolds number. In 
many cases a random perturbation was also added in the 
initial conditions. A run was considered well-resolved if 
the maximum of the current density spectrum was suffi- 
ciently larger than the value of the current density spec- 
trum at the smallest resolved scale (typically by an order 
of magnitude). 

The strongest limitation in these runs turned out to 
be the time of integration. Although large magnetic 
Reynolds numbers can be be achieved in this limit, with 
growth rates of the order of the inverse turnover time, 
the magnetic field with large Prandtl number has much 
longer memory time than a typical turbulent flow. Thus 
to achieve saturation for saturation each run has to be 
integrated for time scales of the order of the diffusion 
time. 



III. THE EXAMINED FLOWS 

Three different forcing functions F were used to gener- 
ate the dynamo flows. The first two are members of the 
ABC family: 



u x = A sm{kfz) + B cos(fcyy), 
u y = C sin(fc^x) + Acos(kfz), 
u z = B sm(kfy) + C cos(kfx). 



(7) 



namo theory 
(V x u = — kju 



The ABC flows have been extensively investigated in the 
literature (l4h 16 , especially for their relevance in dy- 
17-|20|. They are fully-helical satisfying 
Here the case for A = B = C = 1/ v3 
is investigated for kf = 1 and for the case kf = 2. 

The third flow under investigation is a variation of the 
ABC family where only the cosine terms are kept: 



Acos(k f y), 
y — B cos(fc/2;), 



(8) 



C cos(kfx). 



This flow was first investigated in [2lj, l22[ and was shown 
to result in dynamo with large saturation amplitudes of 
the magnetic field p3l - l27j . This flow will be referred as 
the Archontis flow. Unlike the ABC flows, it has zero 
global helicity. In the performed runs in this work have 
A = B = C = v/2/3 and kf = 1. 
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FIG. 1: Time averaged magnetic (squares) and kinetic (dia- 
monds) energy as a function of R m for the Archontis flow. 



IV. RESULTS 



Archontis Flow 



We begin with the non-helical Archontis flow. Fig- 
ure [1] shows the saturation values of the magnetic en- 
ergy Eb = 5l|b|| 2 (squares) and of the kinetic energy 
E u = ^||u|| 2 (diamonds) as a function of the magnetic 
Reynolds number R m . For small values of R m there is 
no dynamo and thus ||b|| 2 = and ||u|| 2 = 1. The onset 
of the dynamo instability is at R mc ~ 42.5. For small 
deviations from this critical value the magnetic energy 
grows as Eb ~ (Rm — Rmc) as a normal mode expansion 
would predict [28|, [2|| , and the kinetic energy sharply de- 
creases. The dynamo is stationary at this stage. As R m 
is further increased the system goes through a series of 
bifurcations that we do not fully explore here. There is 
an other critical value R' mc ~ 80 after which the evolution 
of the magnetic energy appears to be chaotic. We avoid 
using the term turbulence since the flow is in the small 
Re regime, however many common features with turbu- 
lence like power-law scalings are present in this regime. 
It is noted that this chaotic state does not result from 
a random forcing or hydrodynamic turbulence but it is 
solely generated by magnetic instabilities. The dynamo 
remains in this chaotic state for larger R m for all the ex- 
amined runs. Magnetic field energy continues to grow as 
R m is further increased. This increase however appears 
to be logarithmic Eb ~ ln(i? m ) or possibly as a weak 
power law E n ~ R^ with 5 < 0.15. On the other hand, 
the kinetic energy is decreased approaching zero as R m 
is increased. 

Figure [5] shows the Ohmic (squares) and the viscous 
(diamonds) dissipation rate as a function of R m . As the 
dynamo threshold is crossed the viscous dissipation rate 
is decreased and the Ohmic dissipation rate is increased. 
The total dissipation rate is smaller than the laminar 
no-dynamo dissipation rate (ie the system is less dis- 
sipative when the dynamo is turned on). The increase 
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FIG. 2: Time averaged Ohmic (squares) and viscous (dia- 
monds) dissipation rates as a function of R m for the Archontis 
flow. 
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FIG. 3: Magnetic energy spectrum for three values of i? m : 
R m = 2000 (dotted), Rm = 5000 (dashed), R m = 10000 
(solid) for the Archontis flow. The lower solid line shows the 
kinetic energy spectrum for R m — 10000. 

of the Ohmic dissipation rate continues until the criti- 
cal value R' mc , where the chaotic behavior starts. After 
this value both the viscous and the Ohmic dissipation 
rate are decreasing and become almost equal. This de- 
crease continued up until the highest examined magnetic 
Reynolds number. Here, the results of the numerical sim- 
ulations, leave open the possibility that at infinite R m a 
zero-injection, zero-dissipation state could be reached. 

Finally, in figure [3] the energy spectra for the three 
largest examined values of R m , arc shown. As expected 
the cut-off wave number due to magnetic diffusivity is in- 
creased as R m is increased. What is also observed is that 
the bulk of the magnetic energy also moves from the large 
to the small scales as the magnetic Reynolds number is 
increased, with the slope of the spectrum changing from 
a negative value to almost flat. It should also be noted 
that the large scale magnetic energy is decreasing as R m 
is increased while the total magnetic energy is increased. 

Figure S] shows magnetic energy density isosurfaces, for 
the Archontis Flow, R m = 10 4 . Most of the energy is con- 
secrated in elongated almost two-dimensional structures. 



FIG. 4: Magnetic energy density isosurfaces, for the Archontis 
Flow, R m = 10 4 . The surfaces correspond to 50% of the 
maximum magnetic energy density or 10 times the average 
magnetic energy. 



These structures have been reported before in kinematic 
dynamo simulations (23j and are referred to as "magnetic 
ribbons" . Although these structures occupy a small frac- 
tion of the total volume they are responsible for almost 
bringing to a halt the flow in the entire domain. 



B. ABC Flow, k f = 1 

Next we examine the helical ABC flow forced at 
wavenumber kf = 1. This is perhaps the most examined 
flow in periodic boxes for dynamo action. The onset of 
instability and the kinematic regim e of the dynamo has 
been examined in the works [17H19I | . It has been shown 
that the first onset of the dynamo is at R m d — 15.5. This 
mode however stops being a dynamo when R m becomes 
larger than R m c2 — 31. A second unstable mode appears 
at R m c3 = 43 and the dynamo instability is present for 
all larger values of R m . 

Figure O shows the saturation levels of the magnetic 
(squares) and kinetic (diamonds) energy as a function of 
the magnetic Reynolds number. The dynamo instability 
begins at R mc \ as has been previously found and it ex- 
hibits a supercritical bifurcation. At the no-dynamo win- 
dow between R mc i and R m cz the instability is sub-critical 
at both ends of the window. This is revealed by the fi- 
nite amplitude of saturation of the magnetic energy right 
at the onset of the instability. The extend of the sub- 
criticality appears however to be very small: Non-linear 
solutions with finite magnetic energy were found inside 
the no-dynamo window but very close to the boundaries 
of the window. 

For larger values of R m the magnetic energy continues 
to grow and this flow also exhibits chaotic behavior. As 
in the Archontis Flow the magnetic energy shows a loga- 
rithmic (or a weak power law) increase with R rn . Unlike 
the Archontis flow however the kinetic energy saturates 
to a finite value E u ~ 0.125. This is mostly due to the 
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FIG. 5: Time averaged magnetic (squares) and kinetic (di- 
amonds) energy as a function of R m for the ABC, kf = 1 
flow. 
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FIG. 7: Magnetic energy spectrum for three values of R m : 
R m = 2000 (dotted), R m = 5000 (dashed), R m = 10000 
(solid) for the ABC, kf — 1 flow. The lower solid line shows 
the kinetic energy spectrum for R m = 10000. 




FIG. 6: Time averaged Ohmic (squares) and viscous (dia- 
monds) dissipation rates as a function of R m for the ABC, 
kf = 1 flow. 

oscillatory nature of the ABC dynamo that is driving a 
time dependent flow in the non-linear regime. 

The Ohmic (squares) and viscous (diamonds) energy 
dissipation for the ABC dynamo is shown in figure [6] as 
a function of R m . After the initial increase of the Ohmic 
dissipation at larger values of R m , the two dissipations 
become almost independent on R m with the viscous dis- 
sipation being twice larger than the Ohmic dissipation. 
As in the previously examined flow, in the presence of 
dynamo the system is less dissipative. 

The magnetic energy spectra shown in figure [7] appear 
to have a positive slope at small wave numbers. This 
is in contrast with the MHD simulations with order one 
Prandtl number for which at the saturated state most 
of the magnetic energy is concentrated close to the forc- 
ing scale |3Cj Hl|. The positive slope extends to larger 
wavenumbers as the magnetic Reynolds number is in- 
creased. 

The magnetic energy is consecrated in elongated struc- 
tures as shown in figure[8] A grey-scale image of magnetic 
energy through the mid plane of the box shown in figjS] 



FIG. 8: Magnetic energy density isosurfaces, for the ABC 
Flow kf = 1, R m = 10 4 . The surfaces correspond to 50% of 
the maximum magnetic energy density or 10 times the average 
magnetic energy. 

indicates more clearly that these structures are also flat 
and part of sheet-like structures that intensify in regions 
of strong shear; much like the the "ribbons" observed in 
the Archontis flow. 



C. ABC flow, k s = 2 

Finally we examined the case for which the ABC flow 
is forced at wave number kf = 2. The difference with 
the kf = 1 case is that more magnetic modes are present 
in the system that move the onset of the instability to 
smaller values, and allow magnetic helicity to concentrate 
in the large scales. 

Figure [10] shows the magnetic (squares) and kinetic 
(diamonds) energy at saturation for this flow. The onset 
of the instability appears at R m ~ 6 and there is no no- 
dynamo window as in the kf = 1 case. For large values 
of R m the magnetic energy can be seen to increase at 
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FIG. 9: Grey-scale image of magnetic energy through the 
mid plane z = tt of figure [8] (ABC Flow kf = 1, R m = 10 4 ). 
Bright regions indicate large magnetic energy. 
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FIG. 10: Time averaged magnetic (squares) and kinetic (di- 
amonds) energy as a function of R m for the ABC, kf — 2 
flow. 



least as fast as the logarithm of the magnetic Reynolds 
number. This increase is more pronounced than the pre- 
viously examined flows. The kinetic energy saturates at 
a finite value roughly equal to one fifth of the value in 
the absence of dynamo. The two dissipation rates Ohmic 
(squares) and viscous (diamonds) shown in figure HT1 ap- 
pear to be proportional for sufficiently large R m with the 
Ohmic dissipation rate being roughly four times smaller 
than the viscous dissipation. 

The most pronounced difference from the previous 
cases can be seen in the magnetic energy spectra (Fig. 
[T2)) . Unlike the previously examined cases the bulk of 
the magnetic energy is concentrated at the large scales 
as a result of the inverse cascade of the magnetic helicity 
[32l - l37| . In the small scales the energy spectrum appears 
to follow a power law behavior Eb(k) ~ k a with an ex- 
ponent close to a = — 1. 

The magnetic and kinetic helicity as a function of R m 
is shown in figure Q2J The last three points in this fig- 
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FIG. 11: Time averaged Ohmic (squares) and viscous (dia- 
monds) dissipation rates as a function of R m for the ABC, 
kf — 2 flow. 
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FIG. 12: Magnetic energy spectrum for three values of R m : 
R m = 2000 (dotted), R m = 5000 (dashed), R m = 10000 
(solid) for the ABC, kf — 2 flow. The lower solid line shows 
the kinetic energy spectrum for R, n = 10000. 



ure were obtained from the numerical simulations after 
roughly one diffusion time where the magnetic helicity 
was still slightly increasing with time. To obtain a well 
converged averaged value many diffusion times would be 
required that is not feasible with the present computa- 
tional power. Thus for these points the saturated value 
of the magnetic helicity is expected to be slightly higher 
than what is shown in figure [T2J From the figure however 
one can see that the saturated value of the magnetic he- 
licity is slowly increasing with R m while hydrodynamic 
helicity is approaching an asymptotic value. As discussed 
in the introduction this positive value of magnetic helicity 
implies that negative helicity has cascaded to the small 
scales where it was dissipated while positive helicity was 
concentrated at the large scales. At the non-linear stage 
the flow is altered so that although still helical (see fig. 
[T"3| it is less efficient at generating large scale helicity. 
The spectrum of the magnetic helicity is shown in 1141 
Unlike the kinematic stage in which negative helicity is 
concentrated in the small scales and positive in the large, 
at saturation helicity has the same sign in both large and 
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FIG. 13: Time averaged magnetic (squares) and kinetic (di- 
amonds) helicity as a function of R m for the ABC, fc/ = 2 
flow. 
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FIG. 14: Magnetic helicity spectrum for three values of R m : 
R m = 2000 (dotted), R m = 5000 (dashed), R m = 10000 
(solid) for the ABC, k f = 2 flow. 




FIG. 15: Magnetic energy density isosurfaces, for the Archon- 
tis Flow, R m = 10 4 . The surfaces correspond to 50% of the 
maximum magnetic energy density or 10 times the average 
magnetic energy. 
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FIG. 16: Magnetic energy dissipation wavenumber k D as a 
function of R m . 



small scales. This implies that at the nonlinear stage 
some of the large scale helicity has been cascaded to the 
small scales. The transfer of the large scale helicity to 
the small scales has been observed in MHD simulations 
where a scale to scale transfer of the magnetic helicity 
was investigated |37[- It is also noted that the helicity 
spectrum is very steep. If the magnetic field was fully 
helical at small scales then it would be expected that the 
helicity spectrum would scale like H m ~ E^/k ~ fc~ 2 , 
but in fig. [14] the slope is much steeper (H m ~ fc -4 ) im- 
plying that both signs of helicity arc present in the small 
scales with the positive sign only slightly dominating. 

The magnetic structures that appear for this flow can 
be seen in figure [TS] where magnetic energy density iso- 
surfaces are shown. Again here energy is concentrated 
in long and flat structures (ribbons). Although small 
scale structures are also present here they appear to be 
more "organized" to form large scales patterns. Note the 
twisting of the "ribbons" that indicates the presence of 
helicity. 



V. SCALINGS AT SATURATION 

The three examined flows have some some common 
features that is worth further exploring. One of the fun- 
damental relations from the kinematic theory that has 
been shown to hold for a variety of flows is that during 
the kinematic stage the stretching term that scales like 
~ UB/£ LI (where U,B,£ U are typical velocity, magnetic 
field, and velocity length scale respectably) is balanced 
by the magnetic diffusion term that scales like ~ rj B j £ 2 D 
(where £ D is the diffusion length scale). This relation 
leads to the prediction that £ D ~ £ u /y / R m . This re- 
lation continues to hold at the saturated stage. Fig [TBI 
shows the dissipation wave number k D ~ 1 / £ D (defined as 
k 2 D = ((Vb) 2 )/(b 2 )) as a function of R m . For sufficiently 
large R m for all flows the dissipation wave number follows 

1 /2 

a scaling close to k D ~ R„{ . (A best fit gives a value of 
the exponent closer to k D ~ i?^ 4 -) The velocity length 
scale on the other hand appears to depend weakly on the 
magnetic Reynolds number. In figure[T7]we plot the ratio 
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FIG. 17: The kinetic energy dissipation wavenumber fc^ as a 
function of R m . 
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FIG. 18: The effective magnetic Reynolds number Rmeff 
(based on the velocity amplitude and lengthscale at the non- 
linear stage) as a function of R m for the three different flows. 



of the velocity wavenumber k 2 = ((Vu) 2 )/(u 2 )) ~ £~ 2 
to the forcing wavenumber kf. k v for the two ABC flows 
appears to be almost independent of R m and close to kf. 
For the Archontis flow the velocity wavenumber is weakly 
increasing with R m . Thus at the nonlinear stage both 
lcngthscales £ D and l v follow a similar scaling with R m 
as they do in the the linear stage. Furthermore this scal- 
ing also explains why the Ohmic dissipation that scales 
like Cb ~ B 2 /R m l 2 D is of the same order with the viscous 
dissipation that scales like e u ~ U 2 /£ 2 . 

Besides the two length scales it is also interesting 
to investigate the effect of saturation on the magnetic 
Reynolds number. An effective magnetic Reynolds num- 
ber at saturation can be defined as R me ff = U n i/k u r), 
where U n i — (2-En) 1 / 2 is the velocity amplitude at the 
nonlinear stage and k n the velocity wave number as de- 
fined in the previous paragraph. If the dynamo sat- 
urates by decreasing the effective Reynolds number at 
the nonlinear stage then it would be expected that as 
R m is increased, Rmeff would approach the onset value 
R m c- In figure [TBI we plot the effective Reynolds num- 
ber as a function of R m . For the ABC flows the rela- 



tion between the two Reynolds numbers is very close to 
linear. For the Archontis flow the Rmeff is increasing 
with R m but weaker than linear and closer to the scaling 
Rmeff ~ Rm - Figure [T51 then implies that all relations 
described in this work for R m hold also for R me ff for 
the ABC flows, and for the Archontis flow provided the 

3/4 

rescaling R me ff ~ Rm is taken into account. Finally 
since R m ef f is increasing with R m saturation is not due 
to a decrease of the amplitude of the velocity. 

Although the magnetic field does not modify signif- 
icantly the amplitude of the flow to saturate the dy- 
namo it does modify its structure. A measure of the 
strength of this modification is given by the amplitude 
of the incompressible projection of the Lorentz force 
Fj X i, = b ■ Vb — VP. Figure shows the amplitude 
H-Pj'x&H = (Fjxb) 1 ^ 2 as a function of Rm. What is ob- 
served is that for moderate values of R m ( 50 < R m < 
1000 ) the amplitude of the Lorentz force is close to the 
forcing amplitude (||F|| = 1). As R m is further increased 
the amplitude of Fj X b seems to increase to larger values 
for the two ABC flows, while a weaker increase is ob- 
served for the Archontis flow. This suggests that at least 
for moderate R m , Fj X b is in balance the body force F. At 
larger R m the effect of the Lorentz force is probably less 
organized due to the presence of fluctuations and larger 
amplitude of H-Fjxbl arc needed to achieve a balance with 
F. 

One then needs to explain how a balance between 
a force generated by a small scale magnetic field (£d) 
and a large scale force (£jj) is achieved. In figure [20] 
we show the relative amplitude of ||-F}xb|| normalized by 
the amplitude of the current density ||j|| = (V x 6 2 ) 1 / 2 
and magnetic field amplitude ||6|| = (b 2 ) 1 / 2 as a func- 
tion of the magnetic Reynolds number. If the two fields 
b and j = Vxb were randomly organized, the quan- 
tity H-Fjx&ll/flb'H \\b\\) should be of order one. If how- 
ever the two fields were aligned or were organized so that 
their cross product is a potential field the relative am- 
plitude of H-Fjx&H can be a lot smaller. In figure [20l the 
relative amplitude of the Lorentz force appears to de- 
crease as a power-law for all flows with exponent close 

— 1/2 

to H-Fjx&H ~ R m ||j|| ||6||. (A best fit gives a value of 
the exponent closer to -0.4). This can be understood in 
the following way. Saturation is expected to be achieved 
when the Lorentz force is of the same order with the 
external forcing F ~ Fj X b = The amplitude of 

Lorentz force is also expected to scale as Fj X b ~ /||i|| ■ \\b\\ 
where / is a prefactor that depends on the alignment 
and structure of the two fields. Since the current den- 
sity amplitude is controlled by the balance of stretching 

and diffusion as seen in figure \TE\ We expect the scaling 
l li 

\j\ ~ IHI-Rm /L. Equating the two relations we obtain 
that for an order one magnetic field the prefactor / must 

— 1/2 

scale like R m ' . This implies that the cross product of 
the two fields j and b must come closer and closer to a 
potential flow as R m is increased. This is achieved by 
forming the flat structures observed in figures 14191 and [15] 
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FIG. 20: The amplitude of the Lorentz force normalized by 
\\j\\ H&ll for the three different flows as a function of R m . 

where both the direction and the gradient of j x b is in 
the direction perpendicular to their surface. Fj X b then 
varies a much longer length scale than that of the mag- 
netic field that is proportional to the curvature of these 
structures and is of the order of the forcing scale. 

The observed deviations from this balance as R m is in- 
creased only reflect that Fj x t ~ F is not satisfied due to 
the presence of fluctuations that require ||-F)x&|| > \\F\\ 
and thus stronger magnetic field to achieve saturation. 
The behavior of these fluctuations due to magnetic in- 
stabilities thus seem to play an important role for the 
saturation at high R rn and need to be studied further. 

VI. SUMMARY AND DISCUSSION 

The nonlinear behavior of the dynamo instability was 
investigated in the infinite Prandtl number limit for three 
different stationary flows. In this limit large values of 
R m can be examined at moderate resolutions and this 
allowed us to investigate how the non-linear behavior of 
the dynamo scales with R m . What was shown for all 
flows is that the magnetic energy is increasing with R m at 



least logarithmically. This increase was most pronounced 
for the ABC flow forced at kf = 2 for which an inverse 
cascade of magnetic helicity was present. The Ohmic 
and viscous dissipation rates on the other hand varied in 
behavior. For the two ABC flows they were approach- 
ing an asymptotic value as R m was increased while for 
the Archontis flow were slowly decreasing R m . For all 
flows however the Ohmic and viscous dissipation rates 
although operating at different scales were of the same 
order as a result of the balance between stretching and 
diffusion. 

The saturation of the dynamo comes from a balance 
between the Lorentz force and the external body force as 
their similar amplitude suggests. At saturation the mag- 
netic field forms elongated flat structures whose thick- 

— 1/2 

ness scales like R m . The typical lengthscale along the 
other directions was of the order of the forcing length- 
scale. These structures were almost force-free with the 
solenoidal projection of the Lorentz force scaling like 
Fjxb ~ ||i|| H&H/V-Rm- This decrease of Fj X b is the re- 
sult of the organization of the magnetic field in these flat 
structures. The cross product of the magnetic field with 
the current in these flat structures results in a field that is 
both pointing and varying in the direction perpendicular 
to the surface, and thus is close to a potential field. In 
that way the magnetic field is successful in balancing the 
external force that differs both in amplitude (compared 
to B 2 /l D ) and length-scale. 

Another common feature that all flows had was that 
the generated structures were unstable to small scale fluc- 
tuations that resulted in a chaotic behavior. Note that 
these fluctuations are due to magnetic instabilities and 
their origin is probably related to reconnection events. 
However, these chaotic fluctuations have not resulted in 
universal spectra. The observed slopes of the spectra 
varied from positive (ABCkf = 1 flow) to almost flat 
spectrum (Archontis flow) to negative slope (ABCkf = 2 
flow) . Thus if these instabilities lead to universal behav- 
ior it must happen for even larger magnetic Reynolds 
numbers. 

Finally, a comment needs to be made for the applica- 
bility of these results to realistic flows. It is the author's 
belief that although the dynamo has been investigated 
in a somehow idealistic limit some insight can be gained. 
In a realistic flow there exist both large sale structured 
flows forced by a thermal or other instability and turbu- 
lent small scale fluctuations as a result of a turbulent cas- 
cade. Although the turbulent fluctuations result in the 
fastest growing modes in a large Prandtl number dynamo 
flow at saturation the magnetic field configuration must 
be such that the Lorentz force prevents further magnetic 
field amplification due to both the turbulent fluctuations, 
and the large forcing scales. The results in this work help 
in understanding the latter behavior. In particular this 
work gives a clear example of how small scale fields can 
organize themselves to come in balance with a large scale 
forcing. 



10 



Zeldovich Ya. B., Ruzmaikin A.A., Sokoloff D.D. "Mag- 
netic Fields in Astrophysics" 1990 Gordon & Breach Sci- 
ence Pub. 

Cattaneo, F., Hughes D.W. & Kim E.J., "Suppression 
of chaos in a simplified nonlinear dynamo model" Phys. 
Rev. Lett. 76, 2057 (1996). 

Maksymczuk J., & Gilbert A.D., "Remarks on the equi- 
libration of high conductivity dynamos" Geophys. Astro- 
phys. Fluid Dynam. 90 127 (1998) 

Cattaneo F. & Tobias, S.M. "Dynamo properties of the 
turbulent velocity field of a saturated dynamo" J. Fluid 
Mech., 621 205 (2009) 

Tilgner, A. & Brandenburg, A. "A growing dynamo from 
a saturated Roberts flow dynamo" Monthly Notices of 
the Royal Astronomical Society, 391 1477 (2008) 
Schrinner, M., Schmitt, D., Cameron, R. & Hoyng, P., 
"Saturation and time dependence of geodynamo models" 
Geophys. J. International, 182 675 (2010) 
Schekochihin, A. A., Cowley, S.C., Taylor, S.F., Ham- 
mett, G.W., Maron, J.L., McWilliams, J.C. "Saturated 
State of the Nonlinear Small-Scale Dynamo" Phys. Rev. 
Lett., 92, 084504 (2004) 

Schekochihin, A. A., Cowley, S.C., Maron, J.L. 



Phys. Rev. Lett. 



, "Self- 
vol. 92 



Similar Turbulent Dynamo' 
064501 (2004) 

Schekochihin, A. A., Maron J. L., Cowley, S.C., 
McWilliams, J.C, "The Small-Scale Structure of Mag- 
netohydrodynamic Turbulence with Large Magnetic 
Prandtl Numbers" Astrophys. J., 576, 806-813. (2002) 
Schekochihin, A. A., Cowley, S.C., Maron, J.L., 
Malyshkin, L., "Structure of small-scale magnetic fields 
in the kinematic dynamo theory" Phys. Rev. E, vol. 65, 
016305 (2002) 

Moffatt H.K., "The degree of knottedness of tangled vor- 
tex lines" J. Fluid Mech. 35, 117 (1969). 
Gomez, D.O., Mininni P.D., and Dmitruk P. "MHD sim- 
ulations and astrophysical applications Original Research 
Article" Adv. Space Res. 35, 899 (2005). 
Gomez, D.O., Mininni P.D., and Dmitruk P. "Paral- 
lel simulations in a turbulent MHD", Physica Scripta, 
T116, 123 (2005). 

Arnold V.I., C. R. Hebd. Seances Acad. Sci. 261, 17 
(1965) 

Dombre, T. Frisch, U., Greene, J.M., H enon, M., Mehr, 
A. & Soward, A.M. "Chaotic streamlines in the ABC 
flows" J. Fluid Mech. 167, 353 (1986) 
Podvigina O.M. & Pouquet A., "On the non-linear sta- 
bility of the 1:1:1 ABC flow" Physica D 75, 471 (1994) 
Arnold V.I. & Korkina E., "The growth of a magnetic eld 
in a three-dimensional steady incompressible flow" Vest. 
Mosk. Un. Ta. Ser. 1, Matem. Mekh., 3 43 (1983) 
Galloway D.J. & U. Frisch U., "Dynamo action in a fam- 
ily of flows with chaotic streamlines" Geophys. Astro- 
phys. Fluid Dynam. 36 53 (1986). 



[19 

[20 
[21 

[22 

[23 

[24 

[25 

[26 

[27 
[28 
[29 

[30 

[31 

[32 

[33 

[34 

[35 

[36 
[37 



Lau Y.T & Finn J.M., "Fast dynamos with nite resistiv- 
ity in steady flows with stagnation points. Phys. Fluids 
B 5 365 (1993). 

Childress S. & Gilbert A.D., "Stretch, Twist, Fold: The 
Fast Dynamo" Springer, New York, (1995) 
Galloway, D.J. & Proctor, M.R.E. "Numerical calcula- 
tions of fast dynamos in smooth velocity fields with real- 
istic diffusion" Nature 356, 691 (1992) 
Archontis, V. "Study of generation and evolution of mag- 
netic elds in stars using 3D MHD simulations of turbulent 
flows" Ph.D. Thesis, Copenhagen University (2000) 
"Numerical simulations of kinematic dynamo action" Ar- 
chontis V., Dorch S.B.F., & Nordlund A., Astron. Astro- 
phys. 397, 393 (2003). 

Dorch, S.B.F. & Archontis, V. "On the saturation of 
astrophysical dynamos: numerical experiments with the 
no-cosines flow". Solar Phys. 224 171 (2004) 
Archontis, V., Dorch, S.B.F. & Nordlund, A. Nonlinear 
MHD dynamo operating at equipartition. Astron. & As- 
trophys. 472, 715 (2007) 

Cameron, R. & Galloway, D. "Saturation properties of 
the Archontis dynamo" M. Not. R. Astr. Soc. 365 735 
(2006) 

Gilbert, A.D., Ponty, Y., Zheligovsky, V. "Dissipative 
structures in a nonlinear dynamo" larXiv : 1 005 . 5259 
Fauve S. & Petrelis F. "Scaling laws of turbulent dy- 
namos" C. R. Physique 8 87 (2007). 
Fauve S. & Petrelis F. "Saturation of the magnetic field 
above the dynamo threshold" Eur. Phys. J. B 22 273 
(2001). 

Mininni, P.D.; Ponty, Y., Montgomery, D.C., Pinton, J- 
F., Politano, H. & Pouquet, A. "Dynamo Regimes with 
a Nonhelical Forcing" Astrophys. J. 626 853 (2005). 
Gomez D.O. & Mininni, P.D. "Direct numerical simula- 
tions of helical dynamo action: MHD and beyond" Non- 
lin. Processes in Geophys. 11 619 (2004). 
Pouquet A., Frisch U.& Leorat J., "Strong MHD turbu- 
lence and the nonlinear dynamo effect". J. Fluid Mech. 
77 321 (1976) 

Gilbert A.D. & Sulem P.L., "On inverse cascades in alpha 
effect dynamos" Geophys. Astrophys. Fluid Dynam. 51 
243 (1990) 

Brandenburg A., "The inverse cascade and nonlinear 
alpha-effect in simulations of isotropic helical hydromag- 
netic turbulence". Astrophys. J. 550, 824 (2001). 
Brandenburg A., Dobler W. & Subramanian K., "Mag- 
netic helicity in stellar dynamos: new numerical experi- 
ments". Astr. Nachr. 323 99 (2002) 

Gilbert A.D., "Magnetic helicity in fast dynamos" Geo- 
phys. Astrophys. Fluid Dynam. 96 135 (2002) 
Alexakis A., Mininni P.D. & Pouquet A. "On the Inverse 
Cascade of Magnetic Helicity", Astrophys. J. 640 335 
(2006) 



